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1. ۴ a paper published under the above title in the Philosophical 
Magazine and Journal of Science for 1892, Lord Kelvin described a 
graphic method for the solution of dynamical problems. We will quote 
from his paper the following account of his method. 

«In dynamical problems regarding the motion of a single particle in 
a plane the following plan is given for drawing any possible path under 
the influence of a force of which the potential is given for every point of 
the plane. Suppose, for example, it is required to find the path of a par- 
ticle projected, with any given velocity, in any given direction through 
any given point P, (see fig. 1)». 

«Calculate the normal component force at this 
point, and divide the square of the velocity by this 
value, to find the radius of curvature of the path at 
that point. Taking this radius on the compasses, 
find the centre of curvature C, in the line P,¿X, per- 
pendicular to the given direction through P,, and de- 
scribe a small arc, P,P,0,, making P,Q, equal to 
about half the length intended for the second arc. 
Calculate the altered velocity for the position Q,, 
according to the potential law, and, as before for P,, 
calculate a fresh radius of curvature for Q, by find- 
ing the normal component force for the altered direc- 
tion of normal, and for the velocity corresponding to 
the position of Q,. With this radius, find the po- 


sition of the centre of curvature, C,, in P,C,Z, the 


Big.) T: 


line of the radius through P,. With this centre of 
` curvature and the fresh radius of curvature, describe an arc P,P,0), 


making POs equal to about half the length intended for the third arc; 


` calculate radius of curvature for position Q,, draw an arc 757,0 Mand 


` continue the procedure». 
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The practical utility of this method is dependent upon the ease or 
difficulty with which the radius of curvature can be calculated. If the 
potential on which the force depends is a somewhat complicated function 
of the codrdinates of the point, the method immediately becomes so 
troublesome to use, that recourse should rather be had to one of the many 
numerical methods of integration for the calculation of the path. 

We intend, however, to demonstrate an extremely simple, approximate 
construction of the radius of curvature by the ad of a number of equipo- 
tential lines, whereby the above difficulty 18 avoided. 


Let the differential equations for the motion of the point be 


SAN zu 
di? 9x 
dy all 
de ay 0 


DENE (ay 
(a) +(e) = 22 
where U is a given function of « and y. Let us further assume that 


we have drawn in advance a number of equipotential lines, 


U=0, U=}, U= 


Let AMB be a part of the path, and M its point of intersection 


with the equipotential line U=2. Further let M N. be the normal, with 


direction towards increasing O, i.e. towards the centre of curvature of 


the path. Jf we designate as D the part of the normal between the two 

nearest equipotential lines, U E and U ee A we have for the 
n n 

radius of curvature ọ in the point M, the approximate formula 


o = pD. 


By the aid of this formula, the radius of curvature may be easily and 
quickly constructed, when the equipotential lines have been drawn once for 
all, The construction of the path of the point by the aid of the radii of 
curvature thereby becomes very simple and easy of accomplishment. We 
will return later to the question of the possibility of modifying the method 


when the radius of curvature is very large, or when the equipotential 
lines lie far from one another. 
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We will first prove the correctness of the approximate formula م‎ = $D. 
If we designate the projection on the normal M N of the force acting on 
the point M, as K,, and the velocity of the point as v, 
know, the mass of the point being equal to 1, 


we have, as we 


ó » oU RE 
According to system I, v? = 2U, and K, = NV where 3N indicates the 


$ 

ur 

2 

E 

5 Fig 2. 

E 

3 derivative of U along the normal M N. This gives 

7-8 

pee Y (II) 
2 e 2U aN 

2 Let us now consider an assumed point M’ on the normal, and call 


its distance from a fixed point on the latter E, § to be reckoned positive 
in the direction MN. If M’ moves, U will be a function of ¥, and if 
this function be developed in series, we easily obtain the desired approxi- 
“mate formula, o = pD. It was in this way, moreover, that I first demon- 


strated this formula. 
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> 
We attain our end most rapidly, however, 1n accordance with Pro- 


i ing & nction of U, and con- 
fessor Runge's suggestion, by regarding 5 as a fu ; 


sidering the corresponding curve, 


E=f(U) 
where U is an abscissa and § an ordinate’. Let A, B and © be the 
— I 7 + T , , 
three points whose abscissee are 7 > 2 and VE and let AA. be 


, 2 7 = 
parallel with the axis of abscissee; then 4A = = and A C =D. 


(AC) 
Fig. 3. 


Thus if, instead of the direction of the tangent in the point B, we 


take as approximate direction AC, we obtain 
n 
tan 8 = — D 
p==D, 


being the angle made by the tangent with the axis of abscisse.‏ م 


On the other hand, tan م‎ = f’ (U), where the value of U is 2 
1 
Now since 
“4 pi = = = Le 
IN 


1 I am indebted to Professor Runge for information respecting this method of procedure. 
At the beginning of March, 1907, I had sent a description of the method with proofs 
to Professor Mittag Leffler, asking if it were known. He sent it to Professor 
Runge in Göttingen, who drew attention to the above simpler proof. 
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; n 
we also obtain for tan ¢ the value 7 If we compare these two expres- 


sions for tan Û, we obtain 


OED.‏ = م 


It is easy to find an expression for the accuracy attained, for, accord- 
. Bing to Taylor’s formula, we have 


A a U م‎ 
ولا‎ + 4U)=f(U)+f"(U) لمعلا ناه‎ (guy? + a Oe 


and 


; E AA f"(U,—6'. AU) 
U — AU) =f(U,) —f' ( U) 4U += 1X2 (40)? = = 1X2 X 3 RE 


| where 6 and § are numbers between o and r. From this we get 


PAU) —F(U,—4U)=2/"(U) 4U+Ef"(U,+0.4U)+f"(U—8'.40] (40). 


4 If we now choose U, == and 4U= = we obtain 
١ Aet 9 1 rr 7 +4 8 rrr p—0' 
ا( ماسم‎ 
A and thus 
e =PD + e, 
where we have put 
1 É 1010 B 
: A Fr EJ) 
Now if f” ZUERA , lies between — m 


and + m, then 


Enea; m. 


N aa N Re A e a 


m | dU dU وى‎ ZU 
‘In order to find m, f” (U) may be expressed by de’ dea FEE 


¡As & = f(U), we have 
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Thus the first three derivatives of the function U along the normal MV 
may be formed, and f”” (U) may then be found by substitution in the 


above formula. 
We see that the error ¢ diminishes very rapidly with increasing 2. 


2. The above method was worked out for the purpose of simpli- 


(ime the study of the integral curves of the system 


d?R E 20 
T * aR 
d’s 20 
ERE La (II) 


dR Y? daz\? 
E 


where Q is the function 


and y is constant. 

One is led to this system, when it is a question of calculating the 
trajectories of electric corpuscles in space under the influence of terrestrial 
magnetism ^. 

One of my assistants, Mr. E. Steenstrup, engineer, has constructed 
under my guidance, after the above method, more than 50 paths of system 
(IN), answering to various values of the constant y. The practical utility 
of the method has thereby been sufficiently tested. The results of these 
constructions will be fully described in Professor Birkeland's report of 
the Norwegian auroral expedition of 1902—3, in the section that treats of 


the mathematical theory of auroras and magnetic perturbations. 


3. We will here briefly relate our experience with regard to the 
employment of our graphic method of finding the paths defined by 
system (I). 


First a number of equipotential curves, 
I 2 
U=o.U=- Uae ee 
n n n 
1 See my paper: „Sur les trajectoires des corpuscules électrisés dans l’espace sous l'action 


du magnétisme terrestre, avec application aux aurores boréales“, Archives des Sciences 
Physiques et Naturelles; July — October, 1907. Geneva, 
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are calculated, and the curves drawn upon paper. n is assumed so large 
that the curves on the diagram are not too far distant from one another; 

> 
in some places it will generally be necessary to assume larger valués for 


n than at other places. For the sake of clearness, it is recommended that 
all curves ur be dotted lines for م‎ odd, and continuous lines for p 


even, and that the odd numbers, 1, 3, 5, 7, etc., be written in the spaces 


between the continuous curves, from U = o onwards, as shown in figure 4. 


Fig. 4. 


If the point M of the path falls upon one of the dotted curves, the 
radius of curvature then equals the portion D of the normal between the 
two nearest continuous curves multiplied by the number of the interval. 
If M falls upon one of the continuous curves, we have to multiply D by 
the even number lying between the numbers of the adjacent intervals. 

The graphic construction of the path may now be accomplished in 
the manner given by Lord Kelvin; but in addition to compasses, it is an 
advantage to use a narrow tape divided into millimetres, by which the 
centre of curvature can be quickly marked. By means of pins, this tape 
is attached to, and detached from, the points in the path and their centres 
of curvature alternately, as the construction progresses. The product 
pD may be found by the aid of a sliding-rule. i 

More accurate results are obtained by substituting transparent gelatine- 
paper, cross-ruled, for the tape. 

Before proceeding to a treatment of the cases in which the method 
must be modified, we will consider the connection between the radius of 
curvature g of an assumed path through the point M (see fig. 5), and the 


radius of curvature e, of a path through the same point, tangent to the 


ae 
y 


y = i 
ue Mr 


